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A.  STATEMENT  OF  PROBLEM 

Predictions  of  rates  and  extent  of  diffusion  by  atmospheric  turbulence  near 
the  ground  surface  have  been  attempted  with  both  Eulerian  and  Lagrangian 
dispersion  models  but  the  success  of  such  models  will  depend  on  how  accurately 
they  mimic  natural  air  flow  patterns.  When  the  land  surface  is  covered  with  tall 
vegetation,  the  air  flow  within  the  canopy  and  immediately  above  takes  on  unique 
characteristics.  These  characteristics  are  initially  revealed  by  the  presence 
of  repeated  patterns  in  time  series  of  scalars  such  as  temperature  or  humidity, 
in  which  the  scalar  slowly  increases  over  a  period  of  a  few  tens  of  seconds  and 
then  abruptly  returns  to  an  earlier  level.  These  patterns  are  usually  referred 
to  as  ramps  and  are  similar  in  appearance  but  not  in  physical  mechanism  to  the 
temperature  patterns  observed  during  the  passage  of  convective  plumes  higher  In 
the  atmospheric  boundary  layer  (Kaimal  and  Businger,  1970;  Wilczak,  1984). 

Examination  of  the  three-dimensional  velocity  field  in  conjunction  with 
individual  ramp  patterns  reveals  a  consistent  and  spatially  coherent  flow 
structure  composed  of  an-  ejection/sweep  combination,  in  which  the  rapid 
temperature  or  humidity  drop  is  nearly  coincident  with  *;he  change  in  direction 
of  the  vertical  velocity  from  updraft  to  downdraft  (Gao  et  al. ,  1989).  •  It  has 
been  shown  that,  for  forests,  coherent  structures  contribute  greatly  to  the  total 
vertical  fluxes  of  heat,  water  vapor,  and  momentum  between  the  vegetation  and 
the  atmosphere.  The  appearance  of  humidity  ramps  in  near-neutral  conditions, 
and  of  Inverted  temperature  ramps  when  the  air  is  statically  stable  makes  it 
very  likely  that  the  structures  are  formed  by  the  high  wind  shear  in  the  reglcn 
near  the  canopy  top,  and  that  they  are  not  only  a  , consequence  of  thcmal 
'convection. 

For  some  time  it  has  been  noted  that  the  structure  of  atmospheric 
turbulence  immediately  above  and  within  vegetation  canopies  has  features  that 
distinguish  it  from  the  flow  observed  over  smooth  terrain.  A  variety  of 
statistical  properties  demonstrate  these  differences,  which  we  belleVe  are 
reflections  of  the  organized  motions  discussed  above.  Spectra  and  cospectra  have 
shapes  that  are  more  sharply  peaked  than  those  over  smooth  terrain;  longitudinal 
and  vertical  velocity  perturbations  are  more  strongly  correlated;  quadrant 
analyses  show  a  dominance  of  sweeps  over  ejections,  and  that  turbulence  Is  highly 
intermittent;  and  longitudinal  and  vertical  velocities  ate  strongly  skewed. 


2 


positively  in  the  former  and  negatively  in  the  latter  case. 

Extensive  data  sets  are  available  to  the  author  from  turbulence  studies 
In  three  distinctly  different  stand  t3rpes:  a  deciduous  forest  from  full  summer 
foliage  to  fully  defoliated  (Shaw  et  al.,  1988),  an  agricultural  field  crop  of 
maize ,  and  two  orchards .  One  of  the  orchards  was  the  almond  orchard  at  Chico , 
California  in  which  part  of  Project  WIND  was  staged.  Temperature  and  humidity 
ramps  had  been  observed  in  all  of  the  data  sets. 

Knowing  that  there  is  a  large  difference  in  the  statistical  properties  nf 
turbulence  between  smooth  and  rough  land  surfaces,  the  quest! jn  arises  whether 
vegetation  of  different  heights  produce  fundamental  differences  in  the  character 
of  the  turbulence,  or  similar  characteristics  but  possibly  scaled  according  to 
the  physical  dimensions  of  the  sy:>tem.  Differences  that  rai^h:  exist  could  appear 
in  both  the  statistical  descriptions  of  the  flow  and  in  the  degree  of 
organization .  of  the  flow  structures.  The  objectives  of  the  study  were  to  inspect 
turbulence  data  from  distinctly  different  stands  of  vegetation,  to  compare  them 
in  terms  of  statistical  properties  of  flow  and  scalar  fields,  and  to  develop  a 
.technique  fcr  examining  organized  motions  and  to  survey  and  compare  the  data 
sets. 


B.  SUMMARY  OF  RESULTS 

The  study  involved  the  analysis  of  four  existing  data  sets.  Three  of  these 
data  sets  were  described  In  the  research  proposal;  namely  a  deciduous  forest 
experiment  held  in  Ontario,  Canada  1986/87  In  which  the  authot*  was  a  participant; 
a  study  of  turbulence  in  an  almond  orchard  near  Chico,  California  in  1987, 
performed  as  part  of  Project  WIND;  and  measurements  in  a  maize  crop  near  the 
University  campus,  Davis  in  1988.  The  fourth  data  sec  arose  from  a  study  held 
in  an  English  walnut  orchard  near  Winters,  California  in  1989  (Zhang  ec  al., 
1991).  The  results  are  presented  , in  . two  sections,  the  first  dealing  with 
statistical  properties  of  the  observed  turbulence  fields,  and  the  second  dealing 
with  conditional  sampling  of  coherent  structures  and  the  appearance  and  detection 
of  scalar  ramps.  All  data  sets  were  broken  Into  30-mindte  blocks  at  the  time 
of  data  logging,  and  statistical  analyses  are  30*mlnuce  averages. 


It  hai  been  known  for  iiome  time  that  the  probability  distribution  of  wind 
velocity  inside  and  immediately  above  plant  canopies  is  much  more  strongly  skewed 
than  the  air  flow  through  the  remainder  of  the  atmospheric?  surface  layer,  or  over 
smooth  surfaces  (Haitani,  1978;  Raupach,  1981;  Shaw  and  Seginer,  1987).  Skewness 
is  a  measure  of  the  asymmetry  of  a  distribution.  It  is  generally  found  that 
longitudinal  velocity  is  positively  skewed  and  vertical  velocity  is  negatively 
skewed  within  the  air  spaces  of  a  rough  surface,  meaning  that  large  downdrafts 
with  large  streamwise  velocity  are  much  more  likely  than  strong  updrafts  of  low 
momentum  fluid.  The  skewnesses  decrease  to  small  values  at  about  twice  the 
height  of  the  vegetation. 

Figure  1  compares  skewness  of  longitudinal  and  vertical  velocity  ca'^culated 
from  data  for  each  of  the  four  locations,  takon  during  weakly  unstable  conditions 
(-0.15<h/L<0,  except  for  the  almond  otchard,  for  which  -0.5<h/L<0,  where  h  is 
stand  height  and  L  is  the  Monin-Obukhov  length).  Skewness  is  a  non-dimensional 
quantity  (the  third  moment  divided  by  the  cube  of  the  standard  deviation)  and 
the  vertical  axis  is  also  normalized  according  to  the  height  of  each  of  the 
canopies.  The  four  results  are  quite  similar  In  form  and  magnitude.  At  about 
twice  canopy  height,  all  skewnesses  are  small.  ■  Vertical  velocity  skewness 
becomes  Increasingly  negative  to  values  of  about  -1  in  the  middle  and  lower  part 
of  the  canopy,  while  streamwise  velocity  skewness  attains  rather  constant, 
positive  values  of  about  0.7  In  the  orchard  and  maize  canopies.  The  forest  is 
slightly  different  in  that  maximiim  longitudinal  velocity  skewness  ^  found  near 
treetop  height  arid  skewness  diminishes  steadily  with  depth  in  the  forest. 
Otherwise,  there, is  quite  remarkable  similarity  between  the  different  sites  when 
the  vertical  dimension  is  normalized  by  the  height  of  each  iruilvldual  canopy. 

While  high  positive  u-skewness  and  negative  w-skewness  Imply  that  strong 
downbursts  of  high  morwiitum  fluid  intermittently  penetrate  the  canopy,  the 
statistic  does  not  establish  that  longitudinal  and  vertical  velocities  are,  in 
fact,  related  to  each  other.  It  is  demonstrated  in  Figure  2,  however,  that  u 
and  w  are  quite  strongly  correlated  in  the  upper  part  of  each  of  the  canopies. 
The  figure  shows>  vertical  profiles  (plotted  against  normalized  height)  of  the 
correlation  coefficients  between  longitudinal  and  vertical  velocity  (r,^)*  and 
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between  vertical  velocity  and  temperature  (r^j)  for  each  of  the  four  sites.  It 
is  apparent  that  the  highest  correlation  coefficients  are  found  in  the  upper 
layers  of  the  canopies  but  the  largest  height  variations  occur  in  the  forest, 
where  a  distinct  peak  is  seen  at  the  first  observation  level  below  treetop 
height.  It  is  possible  that  lack  of  resolution  in  the  instrument  arrays 
prevented  such  a  phenomenon  from  being  observed  in  the  shorter  stands.  In  any 
case,  the  correlation  coefficients  are  considerably  higher  than  those  reported 
over  smooth  terrain.  For  exainple,  Haugen  et  al.  (1971)  reported  r^^  equal  to 
about  -0,3  and  r„T  equal  to  about  0.35,  for  the  same  stability  range  as  our 
study,  in  measurements  over  wheat  stubble  in  Kansas.  Here,  we  find  correlation 
coefficients  of  about  -0.5  and  0.5,  respectively.  These  large  degrees  of ■ 
correlation,  especially  those  between  the  velocity  components,  are  indicative 
of  a  significant  amount  of  organization  in  the  turbulent  flow. 

The  statistical  technique  of  quadrant  analysis  (tflllmarth  and  Lu,  1974; 
Finnigan,  1979;  Raupach,  1981;  Shaw  et  al.,  1983),  was  applied  to  conditionally 
sample  turbulent  longitudinal  and  vertical  velocities  into  quadrants  according 
to  the  combination  of  signs  of  the  two  signals  and  their  contributions  to  the 
'Reynolds  stress  (sweep,:  u'>0,  w'<0;  ejection:  u'<0,  w'>0;  outward  interaction: 
u'>0,  w'>0;  inward  interaction:  u'<0,  w'<0) .  Sweeps  and  ejections  both 
contribute  to  the  downward  flux  of  momentum,  while  the  Interaction  terms  act  in 
the  opposite  sense.  We  found  only  minor  differences  between  the  stands.  All 
showed  the  previously  reported  pattern  in  which  sweeps  exceed  ejections  in  terms 
of  their  contributions  to  the  momentum  flux  in  the  upper,  levels  of  the  canopies, 
and  with  the  interaction  terms  being  small  by  comparison.  While  canopy  density 
appears  to  have  little  Influence  on  velocity  skewnesses,  stress  fractions 
obtained  using  quadrant  a:ialysefi  show  that  the  relative  contributions  by  sweeps 
and  ejections  decrease  with  depth  in  the  forest  mosv  quickly  when  the  fotest  is 
fully  leafed. 

Attempts  to  compare  power  spectra  of  velocities  observed  in  the  different 
stands,  and  to  compare  spectra  with  those  reported  elsewhere  (e .g. ,.  Baldocchi 
and  Meyers,  1988)  have  not  been  very  successful  because  of  differences  and 
difficulties  in  i.ormallzing  spectra.  It  is  normal  to  want  to  eliminate  the 
effects  of  observation  height  and  wind  speed  from  calculated  spectra.  The  most 
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common  method  to  accomplish  this  is  to  form  a  normalized,  and  non-dimensional 
frei^uency  f^nz/U,  where  n  is  the  natural  frequency,  z  Is  the  measurement  level, 
and  U  is  the  mean  wind  speed  at  height  z.  Divi.lcn  by  the  local  wind  speed 
implies  an  acceptance  of  Taylor's  hypothesis  such  that  the  physical  dimension 
of  eddies  can  be  obtained  from  a  time  trace  of  fluctuations  by  assuming  that 
eddies  translate  along  at  a  speed  proportional  to  the  mean  wind  at  that  level. 
If  canopy  turbulence  is  strongly  influenced  by  coherent  structures,  which 
translate  as  entities  with  the  same  velocity  at  all  levels,  the  normalization 
described  above  would  be  inappropriate,  in  such  a  situation,  division  by  the 
structure  translation  velocity,  or  by  wind  speed  at  a  fixed  level,  would  be 
better.  In  at  least  one  case,  we  have  seen  th^t  spectra  at  different  levels 
within  the  canopy  are  best  normalized  by  a  single  wind  speed  and  by  a  single 
height,  because  tne  natural  frequencies  of  the  peaks  of  the  spectra  match  at  all 
obs'-.rvation  levels.  It  is  apparent,  however,  that  if  contributions  to  turbulence 
arise  from  both  vertically  coherent  eddies  and  from  local  processes,  there  is 
no  consistent  method  for  matching  spectra  from  different  heights.  - 

Overall,  our  analyses  show  that  normalized  statistics  are  remarkably 
similar  for  all  four  of  the  stands  when  vertically  scaled  by  the  height  of  the 
top  of  the  vegetation.  This  similarity  exists  despite  the  large  differences  in 
the  dimensions  and  vegetative  structures  of  the  stands. 

Conditional  sampling  analyses 

The  primary  objective  of  this  part  of  the  study  was  to  optimize  a  scheme 
for  detecting  scalar  raepp  patterns  observed  in  traces' of  temperature  or  humidity, 
and  to  survey  the  four  data  sets  for  coherent  structures  so  identified.  Included 
in  this  plan  were  the  goals  to  find  a  way  to  scale  the  frequency  of  appearance 
of  coherent  structures  and  to  characterise.  In  some  quantitative  manner,  ramp 
patterns  in  the  different  data  sets  in  terms  of  their  duration  and  degree  of 
organization;.  We  have  not  yet  completed  this  task.  Ve  have  performed  surveys 
of  Che  data  sets,  as  we  report  below,  but. show  only  early  results  of  our  attempts 
to  apply  wavelet  transforms, to  the  data. 

Figures  3a* d  show  examples  of  traces  bf  temperature  from  each  of  the  four 


6 


studies.  It  is  apparent  that  the  most  clearly  defined  patterns  are  from  the 
talleiit  canopy,  and  we  have  devoted  more  time  to  their  examination  than  to  the 
other  data  sets.  An  algorithm  (Shaw  et  al.  ,  1989)  to  detect  the  rapid  scalar 
change  at  the  termination  of  each  ramp  by  (i)  moving  through  the  data  to  seek 
a  relatively  large  change  in  the  scalar  occurring  over  a  short  interval  of  time 
(of  the  order  of  3s  for  the  forest  data),  ."nd  (ii)  Immediately  examining  the 
scalar  at  a  series  of  lover  observation  levels  for  a  continuation  of  the  change, 
proved  to  be  successful.  Applied  to  the  forest  data,  this  multilevel  scheme 
showed  that  a  reasonably  consistent  number  of  ramps  (of  the  order  of  1  per  100s) 
were  detected  in  the  htimidity  traces  during  each  half  hour  run  (Figure  4) , 
despite  the  fact  that  conditions  changed  during  the  observation  period  fiom 
unstable  to  stable.  On  the  other  hand,  temperature  ramps  during  the  same  period 
matched  the  number  of  humidity  ramps  at  the  start  and  end  of  the  period  but,  when 
conditions  were  close  to  thermally  neutral,  few  temperature  ramps  were  detected. 
We  believe  that  stability  effects  on  the  appearance  of  coherent  structures  are 
secondary  but  that,  in  order  for  scalar  ramps  to  appear,  it  Is  necessary  that 
there  be  a  vertical  distribution  of  the  scalar.  During  near-neutral  stability, 
it  is  likely  that  (potential)  temperature  Is  uniformly  distributed  with  height, 
leaving  none  of  the  characteristic  traces  th/,t  remain  on  the  huiildity  traces 
during  a  continuing  transpiration  process.  ' 

The  same  objective  analysis  scheme  was  used  to  identify  ramp  patterns  in 
time  traces  of  temperature  from  all  four  data  sets.  It  became  apparent  that  the 
frequency  of  occurrence  of  ramps  increased  as  canopy  height, decreased  but  that, 
when  frequency  of  occurrence  was  normalized  with  canopy  height  and  wind  speed 
(£h/U),  ramp  frequency  was  greatest  for  the  forest  and  least  for  the  much  shorter 
agricultural  crop.  The  normalized  frequency  decreased  with  increasing  wind  speed 
in  all  eases. 

Clearly,  the  normalization  described  above  Is  inadequate  because  it  does 
not  collapse  the  data  onto  a  single  curve.  A  much  better  relationship  was  found 
When  the  frequency  of  occurrence  was  plotted  against  a  measure  of  Che  vl.^id  shear,, 
such  as  the  ratio  of  the  friction  velocity  to  the  canopy  height  (lu/h).  Figure 
5  shows  data  from  all  four  canopies  and  demonstrates  a  much  tighter  relationship 
when  presented  in  this  fashion.  Notice  that  freqiiency  and  wind  shear  have  the 
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same  dimensions,  suggesting  that  such  a  relationship  is  to  be  expected.  However, 
data  from  the  three  canopy  types:  forest,  orchard,  and  field  crop,  barely 
overlap,  and  data  from  the  forest  alone  appear  to  have  a  steeper  slope  than  those 
from  the  other  canopies.  The  arguments  presented  here  are  weaker  than  they  would 
have  been  had  the  data  from  the  different  sites  interspersed  on  the  graph. 

It  is  evident  that  the  scalar  signals  are  somewhat  "noisier*  in  the  shorter 
canopies  (Fig.  3),  making  ramp  identification  more  uncertain.  It  is  possible 
that  the  sJ.rpe  of  the  curve  in  Fig.  5  decreases  with  increasing  values  of  u*/h 
because,  in  the  shorter  canopies,  the  signal- to-noise  ratio  is  smaller.  This 
might  also  be  the  reason  that  the  cur^e  does  not  apoear  to  extrapolate  to  the 
origin  of  the  graph.  We  found  it  necessary  to  develop  an  objective  scheme  for 
evaluating  scalar  traces  to  reveal  not  only  ramp  frequency  but  also  to  provide 
measures  of  the  duration  of  coherent  structures  and  of  the  signal -to -noise  ratio. 
To  this  end,  we  are  examining  the  wavelet  transform  technique  but,  as  stated 
earlier,  this  part  of  our  study  is  not  yet  complete. 

In  the  wavelet  transform  method  (Daubechies,  1990;  Carnage  and  Mahrt,  1990; 
Liandrat  and  Moret-Bailly,  1990),  a  short  term  correlation  is  formed  between  the 
scalar  signal,  and  a  basis  function,  such  as  the  Haar  function  (for  -T/2<t<0,  the 
function  equals  +1;  for  0>t>T/2,  the  function  exiuals  *1).  The  correlation 
between  the  signal  and  the  basis  function  is  obt  lined  as  a  function  of  time 
progressing  through  ^he  data  record.  An  example  of  a  tem^>eracure  trae^  and  the 
corresponding  correlation  function  is  shown  in  Fig.  6.  Peaks  In  the  correlation 
align  themselves  with  the  termination  of  the  ramps 

The  Hear  function  proved  to  be  successful  in 
since  its  nhape  does  not  Include  the  rising  portion 
sensitive  to  ramp  size  nor  can  it  be  used  to  evaluate  the  contribution  of  ramps 
to  the  total  signal  variance.  We  have  chosen  to  examine  functions  that  are  not 
true  wavelets  but  which  match  fsore  closely  the  oatuerns  seen  In  scalar  traces. 
An  example  of  such  a  false  wavelet  is  a  function. th^ t  consl.'tts  of  a  slow,  steady 
rise  followed  by  rapid  decrease  to  the  initial  v^lut.  Ve  have  applied  this 
transform  to  artificially  created  data  with  prescribed  durations  ef  ramps  and 
prescribed  intervals  between  ramps,  and  to  which  candom  noise  has  been  added* 


detecting  scalar  ramps  but, 
of  a  ramp  pattern,  it  is  not 


The  technique  has  proven  to  be  successful  in  retrieving  the  duration  and  interval 
between  ramps.' 
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Figure  1. 
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Skewness  of  longicudlnal  and  vertical  velocities  wlchln  and  above 
four  different  stands  of  vegetation. 


HEIGHT  Z/H 


Figure  2. 
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Correlation  coefficients  for  u,w  and  w,T  within  and  above-  four 
different  stands  of  vegetation.  The  forest  data  are  from  twb  periods 
with  LAI-l.S  (diamonds)  and  LAI-0.3  (squares). 
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riguro  3«,  traces  of  temperature  fluctuation  (®C)  from  lx  levels  of  observation 
,  within  and  above  an  IS  n  call  deciduous  forest;  (1)  m;  (11)  z«*34 
m;  (111)  z-18  m;  (Iv)  z-15  m;  (v)  z-ll  m;  (vi)  z-6  m. 


R«.  U.  Trace,  of  temperature  fluctuatlou  'CC)  from  three  level,  vlthln  and 
.bove  a  2.6  ,  tall  field  of  malre  (1)  *^.9  »;  (n,  *-2.5 
2-1.7  n.  ^ 
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Figure  4.  Number  of  temperature  (lower  panel)  and  water  vapor  (upper  panel) 
microfronts  per  second  In  a  series  of  lO^mlhute  runs,  using  an 
objective  detection  sicherae.  Heat  and  vapor  flux  measurements 
(positive  upwards)  are  shown  by  the  solid  lines.  The  stability 
parameter  h/L  1«  indicated  by  the  dashed  line,, 
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f  temperature  trace  from  the 
Ing  correlation,  function. 


